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ABSTRACT 

Context. Magnetic fields are good tracers of gas compression by shock waves in the interstellar medium. These can be caused by the 
interaction of star-formation driven outflows from individual star formation sites as described in the chimney model. Integration along 
the line-of-sight and cosmic-ray diffusion may hamper detection of compressed magnetic fields in many cases. 
Aims. We study the magnetic field structure in the central part of the nuclear starburst galaxy NGC253 with spatial resolutions 
between 40 and 150 pc to detect any filamentary emission associated with the nuclear outflow. As the nuclear region is much brighter 
than the rest of the disc we can distinguish this emission from that of the disc. 

Methods. We used radio polarimetric observations with the VLA. New observations at A3 cm with 7'.'5 resolution were combined 
with archive data at AA 20 and 6 cm. We created a map of the rotation measure distribution between AA 6 and 3 cm and compared it 
with a synthetic polarization map. 

Results. We find filamentary radio continuum emission in a geometrical distribution, which we interpret as the boundary of the 
NW nuclear outflow cone seen in projection. The scaleheight of the continuum emission is 150 ± 20 pc, regardless of the observing 
frequency. The equipartition magnetic field strength is 46 ± 10 fiG for the total field and 21 ±5 yuG for the regular field in the filaments. 
We find that the ordered magnetic field is aligned along the filaments, in agreement with amplification due to compression. The 
perpendicular diffusion coefficient across the filaments is k ± = 1.5 x 10 28 cm 2 s' 1 • £(GeV) 5±0 7 . In the SE part of the nuclear outflow 
cone the magnetic field is pointing away from the disc in form of a helix, with an azimuthal component increasing up to at least 
1200 pc height, where it is about equal to the total component. The ordered magnetic field in the disc is anisotropic within a radius of 
2.2 kpc. At larger radii, the large-scale field is regular and of even parity. 

Conclusions. The magnetic filaments indicate an interaction of the nuclear outflow with the interstellar medium. The magnetic field 
is able to collimate the outflow, which can explain the observed small opening angle of as 26°. Owing to the conservation of angular 
momentum by the plasma in the nuclear outflow, the field lines are frozen into the plasma, and they wind up into a helix. Strong 
adiabatic losses of the cosmic-ray electrons in the accelerated outflow can partly explain why the radio luminosity of the nucleus lies 
below the radio-FIR correlation. 

Key words, galaxies: individual: NGC 253 - magnetic fields - methods: observational - ISM: cosmic rays - galaxies: magnetic fields 
- galaxies: ISM 

1 . Introduction but can also influence the dynamics if its energy density is high 

enough. 

Observations of nearby galaxies show that the large-scale 
H | Magnetic fields are an important ingredient of the interstellar magnetic field is a very sensitive tracer for interaction in 
. . i medium. Not only do they facilitate star formation, but they the interstellar medium that is not seen at any other wave- 
also make an important contribution to the energy budget. The length. An example for this are the colliding antenna galax- 
global magnetic field structure of galaxies can be described by ies that show a ma ximum of polarized intensity outside of 
a spiral field in the disc and a poloid al field in the halo, which their optical extents dChvzv & Bec k 2004). Several galaxies in 
takes an X-shape when seen edge-on {Krause 2008; Beckll2009t the Virgo cluster show signs of interaction possibly due to 
iBraun et"aT]l2010t) . Such X-shaped structures may be formed by ram pressure compression as sharp edges in the total power 
the h ydrodynamic Al outflows driv en by star formation in the emission and enhanced polarization d Vollmer et al.ll20O7T, |20 1 0): 
disc foalla Vecchia & Schavell2008l) . On the other hand, a con- IWezgowiec et al.ll2007l) . MHD simulations bv lde Aville3(12000l) 
curring theory is that these magnetic field structures are part of the interstellar medium reveal a highly filamentary disc-halo 
of a large-scale field, possibly generated by a galactic dynamo interface with transient vertical filaments in the gas and magnetic 
dGissinger et al.l20"09tlHanasz et al.ll2009l) . Therefore, it is desir- fields. The filaments could be regarded as the walls of so-called 
able to understand how the dynamics of the interstellar medium chimneys through which hot gas can break out from the disc 
influence the structure of the magnetic field that is frozen into into the halo whe n a sufficient number of supernovae occur in 
the gas. The magnetic field may only be a tracer for interactions a OB association dNorman &Ikeuchill989H M ac Low & Ferraral 



o 

> 

O 

(N 

o 

u 

6 

03 



(N 
> 

in 
in 

(N 

q 

On 
O 



- 1—1 

X 



2 V. Heesen et al.: Helical magnetic fields in the nuclear outflow of NGC 253 

Table 1. VLA observations presented in this paper. 



A 


V 


Array 


Resolution 


TP noise PI noise 


Dynamic 


Date 


Code 


References 


[cm] 


[GHz] 






[/iJybeam 1 


] 


range 








20 


1.49 


A 


1'.'3 x 2'.' 2 


50 




3600 


Jul. 1987 


AU30 


Ulvestad & Antonucci (1997) 


20 


1.56 


B+C 


7'.' 5 


400 




2500 


Sep. 1990 


AC278 


Carilli et al. (1992) 


6 


4.71 


C 


7'.' 5 


100 


40 


8300 


Jan. 1991 


AC278 


this paper 


3* 


8.46 


DnC 


7'.' 5 


50 


16 


13000 


Oct. 2009 


AH899 


this paper 



'The observations at AA 6 and 3 cm were merged with Effelsberg single-dish data. 



119991) . Farther up in the halo the gas flow may be driven by 
the cosmic-ray pressure, which has a larger s caleheight and 
is not so much affected by radiation losse s dlpavichl 119751: 
iBreitschwerdt et al.lll99ltlBreitschwerdtll2008l) . Observations of 
this phenomenon will be left to the upcoming radio telescopes 
that can observe at the very lowest frequencies and are sensitive 
to the aged cosmic-ray electrons. Magnetic fields may, however, 
already be important for outflows in the disc. 

Searches for small-scale magnetic field structures in galaxies 
that trace the filamentary disc-halo interface are rare. Normally 
the radio continuum emission is smooth and can be described 
as two-component expone ntial haloes with a thin ( 300 pc) and a 
thick radio (1.8 kpc) disc (iDumke & Krausdll998l) . Vertical fil- 
aments are observed in M82 with prominent gaps in the con- 
tinuum emission, possibly created b y vertical magnetic fields , 
which prevent horizontal diffusion dReuter et alJ 1 1992L Il994t) . 
But it was not possible to rel ate these minima to a particular out- 
flow cone d Wills et al.ll 19991) . In NGC 5775 vertical filaments in 
the radio spectral index suggest the outflow of young cosmic-ray 
electrons, but they do not coin cide with the observed radio spurs 
with vertical magnetic fields dDuric etaL I fl^lTmimann et al. 
120001) . In NGC 1569 a huge Ha bubble c an be identified as the 
edge of a radio spur dKeplev et alj|2010t) . In NGC 253 the halo 
magnetic field seems to line up tangentially wit h the edge of 
the sup erbubble filled with soft X-ray emission dHeesen et alJ 
I2009bl hereafter Paper II). 

NGC 253 is an archetypical starb urst galaxy that has 
both a superwind and a disc wind dHeckman et all l2000t 
IStrickland etal] 120021: iHeesen et aill2009ah . The superwind is 
traced by the hot X-ray emission seen in the halo, indicat- 
ing the outflow of hot gas, possibly supernova-heated. The 
hot gas is seen in the centre as a pl ume of 600 pc length 
emerging to the s outheastern (SE) side (IStrickland et al.l [2000t 
iBauer et all 120071) . The northwestern (NW) side is on the far 
side of the disc an d thus heavily abso rbed in X-ray emission 
dPietsch et al]|2000l) . IBauer et al.1 d2007l) detected with sensitive 
XMM-Newton observations the X-ray emission of the actual 
fluid in the outflow, whereas IStrickland et alJ d2000l) used the 
superior resolution of the Chandra telescope to resolve the fil- 
aments consisting of cooled material. The speed of the hot gas 
co uld be measured from H g emission lines to be 390kms~ 1 
by ISchulz & Wegner] d!992l) . Recent, optical integral-field unit, 
spectroscopy observations could measure the increase in the 
outflow velocity with inc reasing distance from the nucleus 
dWestmoquette et al.ll2f3Tlh . 

The disc wind that is traced by a radio continuum h alo is 
different from the nuclear outflow. In Hees en et a i] d2009al here- 
after Paper I), we have shown that the cosmic -ray electrons are 
transported with a bulk velocity of about 300 km s _1 in a mainly 
vertical direction over the full extent of the disc. The outflow 
speed varies only slightly as a function of the galactocentric ra- 
dius, further supporting the idea of a disc wind, because it seems 



unrelated to the nuclea r outflow. A high-resolution study of the 
continuum emission bv lZirakashvili & VoUd d2006l) of the inner 
disc showed a very similar outflow speed. 

In Paper II we investigated the magnetic field structure and 
found that it is an axisymmetric spiral in the disc and X-shape d 
in the halo, as in several other edge-on galaxies dKrausell2008h . 
The origin of the X-shapes is still not clear, but in NGC 253 it 
may be related to the interaction of t he disc wind with the super- 
wind. Simulations of a superwind bv lSuchkov et al.l d 19941) show 
the formation of huge bubbles. The outer lobes of the bubbles 
show an X-shaped structure in projection. If the large-scale mag- 
netic field is compressed in these lobes and aligned with them, 
they would also form an X-shaped structure. 

So far, no high-resolution magnetic field maps have been 
available for NGC 253 that would allow us to trace the struc- 
ture of the X-shaped halo field back into the inner disc and 
study its interaction with the nuclear outflow. The magnetic field 
maps shown in Paper II have a resolution of 30". In this paper, 
we present new observations with less than 10" resolution to 
study the magnetic field structure in the inner disc of NGC 253. 
In Sect. |2] we present our observations and data reduction. In 
Sect.|3]we discuss the structure of the radio continuum emission. 
Section [4] contains the analysis of the magnetic field structure 
measured with polarimetry. Section [5] contains the discussion 
and Sect.[6] our conclusions. In this paper we use the distance 
from lKarachentsev et alJ (l2003h of 3.94 Mpc d" = 19.1 pc), a 
positi on angle of 52°, and an inclination angle of 78° dPencd 
119811) . 

2. Observations and data reduction 

2.1. New VLA observations 

We used radio continuum polarimetric observations with the 
Very Large Array (VLA) at A3 cm with a bandwidth of 2 x 
5OMH2B The observations were carried out in October 2009 
on three days amounting to 1 1 h of total observing time. We 
used a hybrid DnC-configuration to get a circular beam at the 
southern declination of NGC 253 (6 = -25°). We selected three 
pointings with about 3.5 h of on-source integration time each. 
The data were ca librated with J0137+331 (3C48) using the 
iBaars et aT, I (fl977l) flux scale, and J0025-260 was used as a sec- 
ondary (phase) calibrator. 

For the data reduction we used the Astronomical Image 
Processing System (AlPSfl For total power (Stokes 7), the dy- 
namic range has to be large because of the 1 Jy nuclear point-like 
source. We used self-calibration first with only phase solutions 
(2x) and subsequently with a solution in phase & amplitude 

1 The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

2 AIPS is free software from NRAO. 
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(2x), which improved our dynamic range. Of the 24 antennas 
in the array, 2 1 were upgraded E VLA antennas and three were 
VLA antennas. Our final map improved significantly when we 
left out the three VLA antennas, which we did not investigate 
any further. The r.m.s. noise level of 50yifJy beairT 1 is a factor 
of 3 higher than the theoretically expected value. We explain 
this with the high dynamic range, where in our map the peak- 
to-noise ratio is 13000. Since the sensitivity is enough for our 
analysis and the polarization is not affected by it, as shown be- 
low, we did not try to improve our map any further. 

The polarization leakage (D-terms) was calibrated with PCAL 
using the secondary calibrator and observing the variation in the 
polarization as a function of the parallactic angle. The source 
intrinsic polarization varies with the parallactic angle, whereas 
the instrumental polarization stays constant, allowing us to sep- 
arate them from each other. The polarization angle is calibrated 
with the known angle of 3C48. The calibrated (u, v)-data were 
inverted and cleaned until we reached a level of 2x the r.m.s. 
noise level. We used a robust weighting (iBriggsll 19951) with ro- 
bust=0, which gave us a beam size of 7'.'5 both in total power and 
in polarizationQThe three pointings were combined with LTESS 
(part of AIPS) where a linear superposition is used with a co rrec- 
tion for the attenuation of the primary beam dBraunlll988h . We 
took the part of the maps into account that was above 7% of the 
normalized amplitude with respect to the centre of the primary 
beam. For the two pointings adjacent to the nucleus, we sub- 
tracted the nucleus in each "snapshot" of short integration time 
to remove instrumental effects as the nucleus moves through the 
primary beam. This procedure gave us an r.m.s. noise level of 
16/uJybeam -1 in Stokes Q and U, in agreement with the the- 
oretically expected value, where the maps are 81argely free of 
spurious emission from the nucleus. 

The (t/,v)-plane is not filled below a baseline of 30 m cor- 
responding to a (u, v)-distance of 1 kA at /13 cm. Therefore, any 
structures with an angular diameter larger than 3' cannot be im- 
aged. To reconstruct the large angular scales, we merged the 
VLA and Effelsberg data using IMERG, which is part of AIPS. 
The Effelsberg maps were presented in the two preceding papers 
of this series. The merging was done by a Fourier transform of 
the two maps to be combined, where an overlap in the (w,v)- 
plane had to be specified as uvrange, so that the observations 
could be matched. We compared the integrated flux densities in 
several separated regions in the Stokes / map both, in the disc 
and halo. The flux densities agreed within 10% in each of the 
regions for the chosen overlap of - 1 2kA. The Q and U maps 
were merged using the same optimum uvrange as found for /. 

2.2. VLA archive data 

Archive observations at A6 cm in the C-configuration are avail- 
able that have not been published so far. They consist of five 
pointings along the major axis with 30min total integration time 
each, but spread over 6 h to have a good coverage in the (u, v)- 
plane. The data were calibrated as described for the A3 cm obser- 
vations in Stokes / and polarization. As we only have one point- 
ing centred on the nucleus there are no sidelobes of spurious 
polarization. We inverted the (m, v)-data using robust weighting 
(robust = 0) and cleaned down to 2x the r.m.s. noise level. The 
map was merged with Effelsberg single-dish data using IMERG, 
with an overlap of 0.2 - 0.8k/l. The r.m.s. noise level in polariza- 
tion is 40 /vJy beam -1 , in agreement with what we would expect. 

3 All resolutions are referred to as the half power beam width 
(HPBW). 



In total power the r.m.s. noise level is 100 yuJy beairT 1 , limited by 
the high dynamic range of 8300 due to the bright nuclear source. 

Observations in the B-configuration at /120 cm were also 
taken from archiv e data. These were already presented in 
ICarilli et al] (1 19921) albeit with lower resolution, because they 
were combined with a more compact configuration. We at- 
tempted a polarization calibration, but this was not possible, in- 
dicated by a large variation in the polarization angle of the pri- 
mary calibrator with baseline and time. For a successful polar- 
ization calibration, the variation should only be a few degrees at 
most. The observations were centred on the nucleus, where one 
pointing is enough to cover the central part of the galaxy. We 
combined the /120 cm observations in B-configuration in total 
intensity with those in C-configuration, so that we are sensitive 
to angular scales of up to 15', i.e. to all but the largest angular 
scales. The r.m.s. noise level is, wi th 400 jjjy beairT 1 , slightly 
higher than in the map presented bv lCarilli et al] (Il992l) . which 
has a noise level of 350 // Jy beam -1 at 30" resolution (Paper I). 

Furthermore, there is some high-resolution data at 
/120 cm in VLA A-config uration that was presented in 
lUlvestad & Antonuccil d 1997b . although they did not study the 
extended emission. The data were calibrated in the standard 
fashion in AIPS where we used a model for the primary 
calibrator 3C48, because it is resolved on the longest baselines. 
We cleaned the (u, v)-data using robust weighting (robust=0) to 
2x the r.m.s. noise level. The final noise level is 50juJy/beam 
at a resolution of 2'.'2 x l'/3 (/?.a.=-2°), which agrees with 
lUlvestad & Antonuccil d!997l) . We attempted a polarization 
calibration of the data, but it failed again, as indicated by the 
large variation in the primary calibrator's polarization angle. 
The details of the maps presented in this paper are summarized 
in Tabled] 

For the analysis of the rotation measure we also created nat- 
ural weighted maps at AA 6 cm and 3 cm with 1075 resolution. 
At Ab cm we combined the data in polarization from our previ- 
ously published observations in D-configuration with those in C- 
configuration. At both wavelengths the maps were merged with 
Effelsberg single-dish data. 

3. Morphology of the radio continuum emission 

3.1. Nuclear region and bar 

In Fig. Q] we show the /120 cm total p ower radio continuum dis- 
tribution overlaid on an Ha image bv lHoopes et al] d!996h . The 
radio continuum (though mostly non-thermal) traces the H II re- 
gions well, which are visible as clumps in the Ha image. They 
form a ring around the centre, where the spiral arm NW of the 
nucleus is prominent. We discover a long structure in the inner 
region, extending both NE and SW of the nucleus, resembling 
a radio continuum bar that has so far not been reported for this 
galaxy. This radio bar has a position angle of 63° + 1°. In Fig. [2] 
we overlay ed the A6 cm continuum emi ssion onto a 2MAS S K- 
band image that shows the stellar bar dJarrett et al. 1 12001 . The 
radio and the stellar bar have about the same length, and they 
align well, but the former is shifted clockwise with respect to 
the latter (the stellar bar has a position angle of 70° + 1°). 

The galaxy rotates clockwise with trailing spiral arms. 
According to HI observations in the barred galaxies NGC 1097 
and NGC 1365 bar pat t erns rotate more slowly than the gas 
(lOndrechen et al.1 Il989t lOndrechen & van der Huistl Il989l) . so 
that the gas flow enters the bar from behind; i.e. the bar is over- 
taken by the gas. Therefore the radio bar is downstream of the 
stellar bar and on the leading edge with respect to the rotation. 
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Fig. 1. Total power radio continuum at A20 cm with 7'.' 5 resolu- 
tion. Contours are at 3, 6, 8, 10, 12, 15, 20, 30, 50, 100, 200, 
400, 800, a nd 1600 x 400 uJy/b eam. The background is the Ha 
image from lHoopes et all dl996l) . 



This is the same behaviour as in the prominent barred galaxy 
NGC 1365 that has be en extensively studied in radio contin- 
uum dBeck et al.l 120051) . Gas a nd dust usually accumulates in 
the do wnstream region of bars dAfhanassoulalll992l) . ISorai et al.l 
(120001) indeed find that the molecular gas in NGC 253 is located 
in the downstream region, at p. a. = 66° + 1°, so that it is lo- 
cated between the stellar and the radio bar. The inner part of the 
CO-bar of abou t 1 ' ( 1 . 1 kpc) length is also visible in the map by 
ISakamoto etal] d2006l) (Fig.|5j. 

The concentration of gas and dust may be due to a shock 
in the bar's potential that also compresses the magnetic field. 
Compression likely plays a role at least in the eastern bar be- 
cause there the polarization degree is 22 ± 7% (at A3 cm), a high 
value for a location near the midplane. The sequence from down- 
stream to upstream is therefore regular magnetic fields, turbulent 
magnetic field, molecular gas, and stars. This is the reversed or- 
der that was found in the arms of spiral galaxies, as in M5 1 by 
iPatrikeeve t al. (2006J), so it warrants further investigation. 

3.2. Filamentary emission 

Other distinctive morphological features are the two horn-like 
structures that are sticking out of the central region roughly per- 
pendicular to the direction of the bar, each consisting of two fil- 
amentfl The northern one can be clearly seen in the A3 cm map 
that we present in Fig. [3] with a length of 30" = 600 pc. Both 
structures are marked in Fig. [5] as Fl and F2 north and as F3 
and F4 south of the galactic plane. They are emerging from the 
central region that contains the nucleus, which is a strong total 
power source. The nuclear region is extended in the direction of 
the bar and has a length of 1200 pc at this resolution. The nuclear 
source itself is not resolved at 7'.' 5 (150pc) angular resolution. In 
the optical regime the morphology is dominated by absorption 
owing to dust filaments as sho wn in the high-resolu tion HST 
image from the ANGST survey dDalcanton et al.120071) (Fig. [6}. 

4 We note that the observed filaments are apparent structures ob- 
served in projection. They will be explained later as the walls of the 
nuclear outflow cones. 
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Fig. 2. Total power radio continuum at A6 cm with 1". 5 resolu- 
tion. Contours are at 5, 7, 10, 15, 20, 30, 50, 100, 200, 400, 800, 
1600, 3200, and 6400 X 100/dy/beam. The vectors show the 
orientation of the magnetic field (not corrected for Faraday rota- 
tion) where a vector length of 1" is equivalent to 12.5 yuJy/beam 
of polarized intensi ty. The backgroun d is a K band image from 
the 2MASS survey djarrett et al.ll2003l) . 
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Fig. 3. Total power radio continuum at A3 cm at 7'.' 5 resolution. 
Contours are at 3, 5, 7, 10, 15, 20, 30, 50, 100, 200, 400, 800, 
1600, 3200, and 6400 x 50yuJy/beam. Vectors show the orien- 
tation of the magnetic field (assuming negligible Faraday rota- 
tion) where 1" vector length is equivalent to 12.5/t/Jy beam -1 
polarized intensity. The background is the Ho- image from 
IWestmoquette et all d201 lb . which was scaled with a square-root 
function to enhance weak emission. 



The horn-like structures are also visible in the /120 cm A- 
configuration map at a high resolution of 2'.'2 X 1'.'3, shown in 
Fig. [4] as an overlay on X-ray emission from Chandra data (M. 
Hardcastle, priv. com.). At this high resolution the nucleus (RA. 
00 h 47 m 33 s .12, Dec. -25°17'17'.'3) is resolved. The nuclear re- 
gion now shows up as a structure of 900 x 600 pc in extent, with 
the longer extension aligned with the molecular bar, rather than 
with the major axis. Noteworthy is the apparent anti-correlation 
between the X-ray and the radio continuum emission in locations 
where the X-ray emission traces an outflow of hot gas. This is 
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RIGHT ASCENSION (J2000) 

Fig. 4. Total power radio continuum at /120 cm with 1'.'3 x 2'.'2 
resolution. Contours are at -6, -3, 3, 6, 12, 25, 50, 100, 200, 
400, 800, 1600, and 3200 x 50 ji/Jy beam 1 . The background is a 
map of the X-ray emission between 0.5 and 5 keV measured with 
the Chandra satellite convolved with a Gaussian to 2" resolution 
(M. Hardcastle, priv. com.). 

apparent in the thin X-ray filament (in green, R.A. 00 h 47 m 33!2, 
Dec. -25° 17' 10") pointing to the N and in the SE outflow cone, 
where the radio contours bend slightly inwards. The decrease 
in scaleheight indicates enhanced synchrotron loss due to the 
strong magnetic field in the outflow cone (see Paper II for de- 
tails). 

The filamentary structures are better visible in the grey- 
scale map presented in Fig. [5] where the contour lines show 
12 CO J— 2-1 emission observed with the Submillimeter Array 
dSakamoto et alJl2006l) . There are four filaments that all have 
widths of about one beam diameter (40 pc), so that their true 
width is smaller and have a length of up to 400 pc. The filaments 
are not exactly perpendicular to the extension but are slightly in- 
clined, forming rather the boundary of a cone-like structure. The 
position angle of the filaments are -53° ±1° (Fl), -26° + 1° (F2), 
141° + 2° (F3), and 116° + 4° (F4). The opening angles of the 
cones (between Fl and F2) are therefore 27° ± 1° in the NW and 
between F3 and F4 25° ± 4° in the SE, respectively. 

3.3. Outflow cones 

The nuclear outflow of the starburst galaxy is visible as plumes 
of Ha and soft X-ray emission. In the NW on the far side of 
the halo, the emission is heavily absorbed, so that the outflow is 
much less visible at these wavelengths. The Ha filaments bor- 
dering the SE outflow cone coincide in shape and position with 
the radio filaments F3 and F4. The NW outflow cone is seen in 
Ha only as a weak triangular shape, where the border again co- 
incides with the radio filaments, here Fl and F2. In Fig. [3] the 
Ha mission can b e understood as limb brigh tening of the cone 
(see also Fig. 2 in IWestmoquette et al.ll2.01 lb . The Ha emission 
is presumably material around the rim of the cone, which cooled 
from the hot gas in the inner part of the cone with a temperature 
of 10 7 K. The filaments Fl and F2 in the NW are not visible in 
Ha, probably due to absorption. The soft X-ray emission also 
shows some filamentary emissio n at the boundary of the bubbles 
interpreted as cooled material dStrickland et al.l 120001) . Using 
more sensitive XMM observations albeit with coarse resolution, 
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Fig. 5. 12 CO J=2-l emission from observations with the 
Submillimeter Array a t 4" 6 x 2"4 resolution (p. a. = -41°2)from 
ISakamoto etai](l2006l) . Contours are at 3, 6, 12, 25, 50, 100, and 
200 x 7 Jy beam -1 km s~' . The background image is the A20 cm 
total power emission with l'/3 x2'/2 resolution as shown in Fig. [4] 
which shows the weak emission. The green lines show the lo- 
cation of the profiles discussed in Sect. 13.41 The four filaments 
F1^4 and the minimum Ml are denoted. 
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Fig. 6. Magnetic field orientation at A3 cm with 7'.' 5 resolu- 
tion overlaid on an HST image from the ANGST survey 
dDalcanton et al]l2007l) . A vector length of 1" is equivalent to 
a polarized intensity of 21.3 /Jy beam -1 . The image shows only 
the centre and the NE part. 

iBauer et all d2007l) report detecting the actually hot plasma of 
the outflow. 

By comparing the radio contour lines in Fig. [4] with the CO 
contour lines in Fig. [5] we can see several similarities. The ra- 
dio filaments agree in position with extensions in CO emission 
except at the weak filament F4 that is not prominent in CO. 
The minimum Ml on the northern side of the radio is also pro- 
nounced in CO. The similarities in the contour lines hold also for 
the shape for part of the nuclear region. The obvious exceptions 
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Fig. 7. Slices of continuum emission where the dotted lines show the measured profiles and the solid lines the fit to the data, a) Slice 
across the filaments parallel to the central region at a distance of 300 pc. b) and c) Slices along the filaments F2 and Fl, respectively, 
d) Slice along the NW outflow cone. For slices (b)-(d) a constant level of background emission estimated from (a) was subtracted. 



are extensions in CO east and southwest of the nuclear region. 
We will attempt a more quantitative comparison between the ra- 
dio and CO emission in Sect. l3.6l 

3.4. Emission profiles 

A useful way to characterize the distribution of the continuum 
emission is to measure scaleheights. We created profiles of the 
radio continuum emission with a fixed resolution of 7'.'5 = 
150 pc for the filaments Fl and F2. The profiles, along with 
exponential fits, are presented in Fig. [7] where we used a slice 
width of 2", which is smaller than the resolution. Also, we used 
a perpendicular profile across the filaments, which shows the fil- 
aments and the depression between them. For the slices along the 
filaments and the cone (Fig.|7J)-d), we subtracted a constant level 
of background emission estimated from the cross slice (Fig.|7^). 
The contrast between the filaments and the depression depends 
on wavelength. The depression is most prominent at A3 cm, with 
the level of emission in the outflow cone decreasing to the sur- 
rounding level. Therefore, at A3 cm we do not detect any con- 
tinuum emission coming from the interior of the outflow cone. 
At A6 cm the intensity inside the outflow cone is not entirely de- 
creasing to the surrounding level and at /120 cm the depression 
is only shallow. See also Figs.Q][3] where the appearance of the 
filaments clearly depends on wavelength. 



We fitted exponential functions to the emission profiles to 
measure the scaleheights. We found a scaleheight of 150 + 20 pc 
in both filaments at all three wavelengths as shown in Figs. [7]) 
and c. This is remarkable because the scaleheights of the con- 
tinuum emission in the halo increase with increasing wave- 
length (Paper II). This is understood by the increase in the life- 
time of cosmic-ray electrons with increasing wavelength. The 
wavelength-independent continuum scaleheight of the filaments 
indicates that we see a decrease in the magnetic field strength, 
rather than that of the cosmic -ray electrons density. 

We fitted Gaussian profiles to the cross slices in Fig.|7h 
to measure the width of the filaments as a function of wave- 
length. Under the assumption that the filaments have a Gaussian 
profile the observed width can be expressed as FWHM^ bs = 
FWHMgj + (150pc) 2 to account for the limited resolution, which 
is comparable to the width of the filaments. We found FWHMs 
of 330 ± 15 pc 0120 cm), 350 + 15 pc (A6 cm), and 250+ 15 pc 
(A3 cm). We note that these values are considerably higher than 
the filament width of < 40 pc that we inferred from the A2Q cm 
high-resolution map (Sect. l3.2T >. This result suggests a weak, ex- 
tended radio emission component associated with the filaments 
that is not detected in the high-resolution map with a reduced 
sensitivity to extended emission. Both the magnetic field and 
cosmic-ray electrons can contribute to this weak component. 
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Fig. 8. Spectra of the radio continuum flux density as measured 
from the fits to the data shown in Fig. [7] 



The frequency dependence of the filament width suggests a 
transport effect of the cosmic rays (see Sect. 15. lb . A wavelength- 
dependent filament width is also found in supernova remnants, 
where the synchrotron emission in the X-ray regime exhibits ex- 
tremely thin filaments (< 0.1 pc) with strong magnetic fields, 
whereas the radio filaments are much broader as the result of 
the longer lifetime an d diffusion of the lower energy electrons 
dRevnolds et al.ll2011l) . 

3.5. Radio spectral index 

The emission profiles with the same scaleheight in the filaments 
imply that the radio spectral index a (defined as S oc v") is con- 
stant with respect to the distance from the disc. To calculate the 
radio spectral index, we took the flux density from the fits shown 
in Fig. [7] extrapolated to a distance of Opc. For both filaments 
we find a = -1.0 + 0.1, with no significant curvature in the 
spectra shown in Fig. [8] In filament F2 the A3 cm flux density is 
low hinting at strong electron losses. In contrast, the A3 cm flux 
density is high in filament Fl where the thermal fraction may 
be greater. Of course, with only three different wavelengths, we 
cannot measure any curvature precisely but can assume a linear 
spectrum. The radio spectra in the filaments are consistent with 
a non-thermal spectrum within the error interval, with a non- 
thermal radio spectral index of a nt = - 1 . In the interior of the 
cone (see Fig.|7^), the radio spectral index is -2.0 + 0.2 between 
AA 20 and 6 cm, indicating strong energy losses of the electrons. 

In Fig.|9]we present the distribution of the radio spectral in- 
dex between AA 20 and 3 cm. We clipped each map below 4x the 
r.m.s. noise level prior to the combination, so that the maximum 
error is +0.3 in the radio spectral index. The spectrum is flat 
in the centre around the nucleus but steepens significantly away 
from it. The two filaments in the NW and the radio bar show 
up as extensions of the flat central region. The flat spectral in- 
dex of the central region also extends into the radio bar on both 
sides of the nuclear region. At the location of the NW outflow 
cone (R.A. 00 h 47 m 32 s , Dec. -25°16'50"), the radio spectral in- 
dex steepens significantly, as expected for strong radiation losses 
of the electrons (Sect. l5.ll l. We checked carefully that this feature 
is not caused by the merging of the interferometric data with the 
single-dish data. It is already apparent in the radio spectral index 
map constructed from VLA data alone. The steep region south 
of the nucleus (R.A. 00 h 47 m 33\ Dec. -25°17'50"), however, is 
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Fig. 9. Radio spectral index between AA 20 and 3 cm at 7'.' 5 res- 
olution as grey-scale. Contours show the continuum emission at 
A3 cm as in Fig. [3] 




12COJ=2-1 [Jybeam" 1 km s" 1 ] 

Fig. 10. Pixel-by-pixel comparison of the A20 cm radio contin- 
uum intensity with the 12 CO J=2-l intensity at 88 pc resolution. 
The line shows the theoretical expectation using the conversion 
into star formation rates as explained in the text. 



less certain. At A3 cm we observe a minimum in the total power 
emission that may be caused by the high dynamic range close 
to the nucleus. Also the position does not not agree with the SE 
outflow cone (see Fig. [3), so that this feature could be an artefact. 



3.6. Radio continuum as star formation tracer 

The radio continuum emission is a star formation tracer because 
it is sensitive to the thermal and the non-thermal emission. UV- 
photons emitted by massive stars ionize the hydrogen that can 
be seen as thermal Ha and free-free radio emission. Supernova 
remnants that mark the end of life of massive stars accelerate 
cosmic-ray electrons that radiate synchrot ron emission, th e non- 
thermal component of the radio emission. ICondonl(ll992h calcu- 
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Fig. 11. Rotation measure distribution between AA 6 and 3 cm 
with a resolution of 1CX'5. The major and minor axes and the 
position of the sectors are shown. 

lated the star formation rate or, for resolved emission, the star 
formation rate density as 

f SFRD 2 ) =Mxl(T'(J»=a (1) 

where /20 cm is the radio continuum intensity at /120 cm. 

Another star formation tracer is the molecular hydrogen 
traced by the 12 CO J=2-l emission. According to iKennicutl 
(1 1998b the star formation rate density relates to the mass sur- 
face density of combined atomic and molecular hydrogen Zh2 
as 

SFRD ) _ 2 5 x 10 -!0 f ^ ) 1A (2) 
Moyr-^pc- 2 / ' \M Q yr" 1 kpc~ 2 / 

We assume that the entire gas in the nuclear region is of molec- 
ular nature and estimate the gas mass thus from the molecu- 
lar hydrogen alone, which we derived from the CO map. We 
ad opted a ratio betwe en 12 CO(2-l) and 12 CO(1-0) of 0.9 found 
by IWeifi et alj d2005l) for the starburst galaxy M 82. Also we 
used the conversion of 5 x 10 19 molecules cirT 2 /(K kms~') to 
convert from CO to H2 co lumn density, which was derived by 
iDownes & Solomonl d 1 998b for ultraluminous far-infrared (FIR) 
galaxies and is widely used for starburst galaxies in the litera- 
ture. 

In Fig. [10] we compare the radio continuum intensity at 
/120cm and the 12 CO J=2-l intensity on a scale of 476 = 88 pc 
in a pixel-by-pixel plot. The line shows the relation, if we use 
the conversion into star-formation rate densities. The scatter in 
the data points is large, with the radio flux densities on aver- 
age a factor of 4 brighter than theoretically expected. Ott et al. 
J2001 find that the integrated star formation rate in the nucleus 
of NGC 253 as derived from Al, 2 cm radio continuum obser- 
vations agrees with that of the FIR, suggesting that the radio 
emission is a reliable tracer for star formation. Recent studies 
established such a c orrelation in spatially resolved observations. 
iDumas et al.l (1201 ll) confirm a tight correlation between radio 
emission and various other star formation tracers down to a scale 
of 240 pc. We do not find such a tight correlation in the nucleus 
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Fig. 12. Distribution of the depolarization between AA 6 and 3 cm 
at a resolution of 15". Contours show the distribution of the po- 
larized emission at /16 cm with contours at 3, 5, 10, and 20 x 
70juJybeam _1 , the r.m.s. noise level. 

of NGC 253. A possible explanation are transport effects of the 
cosmic rays in the powerful nuclear outflow, which we discuss 
in Sect. 15. 1145.31 If cosmic rays diffuse away from the star for- 
mation sites, the radio emission gets decoupled from other star 
formation tracers, and the relation between star formation rate 
and radio emission breaks down. 

4. Magnetic field structure 

4.1. Magnetic field orientation 

In Figs.|3]and|6]we show the magnetic field orientation at A3 cm, 
where the length of the vectors is proportional to the polarized 
intensity. The maximum polarized intensity is near the nucleus, 
in the adjacent NE region. The distribution is asymmetric with 
respect to the minor axis, the NE half having a higher fraction of 
polarized intensity. There is a local minimum in the NW outflow 
cone where we do not detect polarized emission, coincident with 
the minimum in total power. The magnetic field orientation is 
parallel to the major axis only NE of the nucleus near the stellar 
bar. In all other regions there is a significant vertical component 
particularly in the two radio filaments bordering the NW outflow 
cone. There, the magnetic field is also remarkably well aligned 
with a polarization degree of 15 + 5%. 

For the study of the Faraday rotation we made a map in po- 
larization at A6cm (see Fig.|2j. Because the integration time is 
short (30min), the map is not as sensitive as the A3 cm map. The 
polarized emission is again asymmetric with respect to the minor 
axis with most of the emission in the NE and is aligned roughly 
along the major axis. The magnetic field orientation is different 
in places from that at A3 cm indicating strong Faraday rotation, 
which we discuss in the next section. 

4.2. Rotation measure distribution 

The rotation measure (RM) is given by the line-of-sight com- 
ponent of the large-scale magnetic field. If the RM is positive 
the magnetic field is pointing towards us and if it is negative it 
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points away from us. In Fig.[TT]we present the RM distribution 
between AA 6 and 3 cm at Iff! 5 resolution. The maps of the polar- 
ization angle were clipped below a polarized intensity of 3x the 
r.m.s. noise level, prior to the combination. The RM ambiguity 
is a +n ■ 1200 rad irT 2 . 

The RM changes on scales of a few beam diameters, corre- 
sponding to about 100 pc. This finding agrees with the observa- 
tions in face- on galaxies, which sh ow similar RM distributions, 
e.g. in M51 dFletcher et alJl201 lh . There the halo RM is also 
expected to contribute to the RM signal as a foreground to the 
bright disc. Such frequent RM reversals can hardly be observed 
in edge-on galaxies due to the effect of averaging along the line- 
of-sight. The RM fluctuations may be caused by field reversals 
on small scales, possibly an "anisotropic" field, which is gener- 
ated from an isotropic turbulent field by compression or sheafl 
The RM is typically between -200 and +200radtrT 2 , far 
away from any RM ambiguity between AA 6 and 3 cm, so the 
orientation of the magnetic field vectors at A3 cm are rotated 
by at most ±15°. The rotation angle is small, so that we do 
not apply any correction to the polarization angles at this wave- 
length. We note that the galactic foreground RM for NGC 253 
Z = 97°, b = -88°) is only a small cont ribution at 10 + 5 rad irT 2 
Noutsos et al j2008tlTavlor et al.l2009l) . which we neglect in the 
further analysis. 

4.3. Faraday depolarization 

The observed degree of polarization is usually lower than the 
theoretically highest value of about 74%. Relevant depolariza- 
tion may be caused by wavelength-independent beam depo- 
larization and by wavelength-dependent depolarization effects, 
the Faraday depolarization. The latter consists of differential 
Faraday rotation (related to the regular magnetic field strength 
and thermal electron density within the emitting source, the 
galaxy), of Faraday dispersion related to the turbulent magnetic 
field strength within the emitting source (internal Faraday dis- 
persion), and of a turbulent magnetic field betwee n the galaxy 
and the observer (e xternal Faraday dispersion) dBurnl 1 1966c 
ISokoloff et al.l IT998h . Usually, the differential Faraday rotation 
is expected to contribute most to the observed depolarization 
within a spiral galaxy, because its values are related to the ob- 
served RMs within this galaxy. 

Observationally, the Faraday depolarization is calculated as 
the ratio between the observed degrees of polarization in maps 
at two different wavelengths smoothed to the same angular reso- 
lution or directly by DP = (PF /PI 2 ) x (v 2 /vi) aa , where a nt is the 
non-thermal spectral index. Using a nt = -1 as measured Paper I, 
we calculated DP between AA 6 and 3 cm at 15" resolution pre- 
sented in Fig.[T2l 

The value of DP is for most parts in the galaxy is rather low 
with values between 0.6 and 1.0 where |RM| reaches values up 
to 250 rad rrT 2 . As a result, the expected depolarization at A6 cm 
due to differential Faraday rotation alone leads to 0.65 < DP < 
1 .0 and indicates that other depolarization effects are negligible. 
The depolarization increases with higher distances N of the ma- 
jor axis up to values of about RM = 350 rad irT 2 where the dif- 



5 The observed B-vectors of linearly polarized emission can trace ei- 
ther regular magnetic fields (i.e. preserving their direction within the 
telescope beam) or anisotropic fields (i.e. with multiple field reversals 
within the beam). Anisotropic fields can be generated from turbulent 
fields by shear or compression. To distinguish between these two com- 
ponents, additional Faraday rotation data is needed. The fields observed 
in polarization are called "ordered" throughout this paper. 



ferential depolarization leads to DP = 0.2 , followed by a stripe 
with no polarized emission at /16 cm parallel to the major axis 
but offset in a NW direction. This stripe has an almost linear, 
well-defined boundary, and is probably caused by even stronger 
Faraday depolarization of the regular and turbulent magnetic 
field in the disc that lies in front. 

More to the NW, we again detect polarized emission at A6 cm 
where DP again lies within 0.6 and 1.0. We conclude that if we 
restrict the study of RM to the near side (SE of the major axis), 
our analysis should hardly be affected by depolarization. 

4.4. Disc magnetic field 

According to the mean-field dynamo theory in thin-disc objects, 
the axisymmetric spiral (ASS) mode (m = 0) with even symme - 
try is the dominant pattern of the regular field dBeck etalll996h . 
Our azimuthal RM analysis in Paper II showed a single peri- 
odical variation, in accordance with an m = mode. Our new 
observations allowed us to test this in more detail and to calcu- 
late refined ASS models. We compared the models with the ob- 
servations using averaging in sectors to plot the azimuthal RM 
variation. The radial interval should avoid the starburst region 
and separate the molecular bar region from the disc. The nuclear 
starburst region extends to a radius of « 300 pc. T he molecu- 
lar ba r extends to a galactocentric radius of 3.5 kpc dSorai et al.l 
120001) . while the bright part extends to only about 1 kpc (Fig. [5). 
We chose 2.2 kpc as a compromise. As the outer radius we take 
4.0 kpc, up to where the signal-to-noise ratio allows us the inte- 
gration. In summary, we used three radius intervals, 0.2-1 .0 kpc, 
1.0 - 2.2 kpc and 2.2 - 4.0 kpc, where we chose an inclination 
angle of 50° and an azimuthal spacing of 20° (the position of the 
sector rings is shown in Fig. [Till. For the noise estimate we used 
the r.m.s. noise values of Stokes Q and U. 

Figure[T3]shows that, in the inner and middle ring, |RM| has 
no maxima around the major axis (azimuthal angle 0° and 180°), 
as would be expected for an axisymmetric m = mode. Instead, 
|RM| is strongly enhanced and changes its sign near the minor 
axis (azimuthal angles 90° and 270°) where the m = mode can- 
not contribute to RM. We conclude that the disc magnetic field 
within * 2.2 kpc does not have a regular axisymmetric struc- 
ture. To explain the lack of RM near the major axis, we propose 
that the ordered field in the inner disc is anisotropic, i.e. it con- 
tributes to the polarized intensity but has frequent reversals, so 
that it does not contribute to RM. In Sect. l4.5l we discuss a model 
of the conical outflow that can explain the RM enhancements 
within 2.2 kpc. 

In the outer ring (2.2 - 4.0 kpc), shown in Fig. [14] the even- 
parity axisymmetric disc magnetic field gives a reasonable fit 
to the data, except around 60° - 80°. This deviation can be at- 
tributed to t h e "rad io spur" that has already been identified by 
ICarilli et al] d 19921) in total power radio continuum and which 
has a strong vertical field component. This finding agrees with 
our previous observations (Paper II), where we detected an 
even-parity axisymmetric disc field at lower resolution. Here we 
used the same model parameters as in Paper II, except that the 
anisotropic component in the even disc model is zero. The model 
is therefore no fit to the data. 

4.5. Magnetic field in the outflow cone 

In this section we test whether the magnetic field in the out- 
flow cone can explain the RM enhancements around the minor 
axis seen in Fig. [T3] No polarized emission from an extended 
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Fig. 13. RMs averaged in azimuthal sectors in the inner disc. The red line shows the model of the helical magnetic field in the 
outflow cone as described in the text (no fit). 



Table 2. Parameters for the synthetic polarization maps (other 
parameters as in Paper II). 



Parameter 


Value 


Notes 


Filament width 


40 pc 


from A20 cm 


Base length 


300 pc 


from A20 cm 


Disc field scaleheight /?g lsc 


1.6 kpc 


thin radio disc 


Disc field strength SjJ isc 


4.4yLiG 


from equipartition 


Filament field strength B^ 1 


20 //G 


from equipartition 


Electron density n e 


2 cm -3 


from Ho- and soft X-ray 


Electron scaleheight h e 


1400 pc 


from Ha 


Cone opening angle 2/3 


45° 




CRE scaleheight disc 


800 pc 


thin radio disc 


CRE scaleheight filament 


200 pc 


filaments 


Filament field scaleheight h% 


600 pc 


from equipartition 



Note: CRE = cosmic-ray electrons 



Radius = 2.2 - 4.0 kpc 
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Fig. 14. Azimuthally averaged RMs in the outer disc where the 
red line shows an even-parity axisymmetric spiral model of the 
disc field. 



halo field is detected with our observations as the magnetic field 
vectors are mostly parallel to the disc. The magnetic field only 
opens up farther away from the galactic midplane and at larger 
galactocentric radii (see Paper II). However, our high-resolution 
image (Pig. [5j suggests that the outflow cone has a strong mag- 
netic field aligned with the walls of the cone, so can influence 
RM, particularly in the SE cone that is located in front of the 
bright polarized disc emission and may act as a Faraday screen. 
We made models using an anisotropic magnetic field in the disc 
and a filamentary magnetic field distributed in the walls of the 
conical outflow. In Table[2]we summarize the parameters of the 
model, and we show a sketch in Fig. [15] We assumed a rather 
large cone opening angle of 45° as the cone should open farther 
away from the disc. 

The model RM distribution is shown in Fig. Q~6] along with 
the magnetic field orientation. The magnetic field orientation is 
influenced little by the magnetic field in the outflow cone be- 
cause its |RM| of < 200 rad trT 2 corresponds to a rotation of less 
than 46°. In the NW, the outflow cone is almost not detected in 
RM as it is located behind the bright polarized disc. The verti- 
cal component of the magnetic field points away from the disc 
in the SE outflow cone (see Paper II). However, a vertical field 
alone cannot give rise to the observed strong change in RM. We 



also need an azimuthal component parallel to the rotation direc- 
tion of the disc. Both components together form a helical mag- 
netic field, which generates the strongest RM signals in the walls 
of the outflow cone. We assumed that the azimuthal component 
rises linearly to a height of 1200pc, where the azimuthal com- 
ponent is equal to the vertical field component (opening angle 
a = 45°, see Appendix). As shown in Fig. [15] the field lines 
are vertical close to the disc but start to wind up with increasing 
height. 

Our model can basically explain the observed RM distribu- 
tion, especially the RM jump between 60° and 80° in the mid- 
dle ring (1.0 - 2.2 kpc). The red lines from the model shown 
in Fig.[l3]have the right properties, but are not good fits. In the 
inner ring, the RM peak is at an too small azimuthal angle, prob- 
ably because the cone is not oriented perfectly perpendicular to 
the disc. The RM amplitudes of the model are too small in the 
middle ring, which indicates that the azimuthal field components 
increases faster than linearly with radius. An improved model 
would need more parameters than those given in Table[2] but 
with the present observations we did not attempt to constrain 
them. 
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Fig. 16. Model RM distribution and magnetic field orientation at 
A3 cm at 1075 resolution as explained in the text. 



Fig. 15. Sketch of the magnetic field in the walls bordering the 
SE outflow cone, which appears as filaments when seen in pro- 
jection. Solid lines show field lines on the front side of the cone, 
whereas dotted lines show them on the rear side. 



4.6. Magnetic field strengths 

We estimated the magnetic field strength from the radio inten- 
sities by assuming equipartition between the energy densities 
of the magnetic field and the total cosmic rays dBeck & Krausd 
120051) . As we are only measuring the electrons we have to make 
an assumption about the proton to electron ratio, K = 100, the 
usually assumed value for galaxies. 

This choice of the iT-factor is motivated by theory and ob- 
servations. Firstly, if the acceleration process produces a power 
law in momentum of the cosmic rays and the initial numbers 
of protons and electrons are equal, then one can infer that for 
relativistic energies the /T-factor depends o nly on the mass ra- 
tio of protons t o electrons, as shown e .g., bv lSchlickeiserl (120021 
Sect. 19.4) and lBeck & Krausel d2005l their appendix). Secondly, 
local measurements in the solar system can directly lead to 
the proton and elec tron spectrum and corroborate K = 100. 
lAdriani et al.l (1201 ll) present an electron spectrum based on the 
satellite experiments PAMELA. From their Fig. 1 one can in- 
fer th e electron flux as .18 s -1 irT 2 GeV -1 (converting their 
units). [Shikaze et al] (120071) present a proton spectrum based on 
the balloon experiment BESS. From their Fig. 7, a proton flux 



of 20 s m GeV is inferred. Thus, the solar system value 
agrees with K = 100. 

The ^-factor is nevertheless somewhat uncertain for exter- 
nal galaxies and thus induces uncertainties in the equipartition 
magnetic field strength. Particularly, K may be higher in case of 
strong energy losses of the electrons, where the electrons have 
died away and only the protons (and nuclei) remain. This may 
be the case for galactic haloes which contain an aged population 
of cosmic -ray electrons. As we have seen in Sect. 13. 51 the syn- 
chrotron losses in the nuclear outflow of NGC 253 are impor- 



tant, so that the /^-factor may well be larger than 100; however 
the magnetic field strength depends on K 1 ^ 4 , so that a 50% error 
in K produces only a 13% error in the magnetic field strength. 
Summarizing, our magnetic field strengths should be taken as 
lower limits. 

The radio intensities were taken at AA 20, 6, and 3 cm 
from the fits to the filaments, as shown in Fig. [7] and from the 
maps in Figs.Q][3] where we subtracted the background emis- 
sion from the disc. The thermal contribution that also needs to 
be subtracted (expect for the nucleus) was assumed to be 10% 
(/120 cm), 20% (/16 cm), and 30% (A3 cm). The nuclear region 
has a size of 600 pc, which we took as the integration length 
along the line-of-sight, resulting in a magnetic field strength 
of 160 + 20 /iG. For the region of the outflow cone we took 
500 pc, which is the diameter of the outflow cone. The result 
is fi to t = 46 ± 10 /iG for the total field and B old = 21 + 5/iG 
for the ordered field. The intensity of the filaments Fl and F2, 
observed at /120 cm at high resolution (Pig.[5j with a width of 
< 40 pc, yields a lower limit of the total field strength of 40 /iG 
and 18 /iG for the ordered field. 



Table 3. Equipartition magnetic field strengths in various re- 
gions. 



Region 




Sold 


Pol. degree 




UiG] 


U>G] 


[%] 


Nuclear region 


160 ± 20 


< 5 


< 0.2 


Filaments F1+F2 


46 ± 10 


21 ±5 


21+7 


Eastern bar 


24 + 5 


11 ±3 


22 ±7 


Western bar 


28 ±6 


5±2 


4± 1 



Note: the polarization degree of the synchrotron emission refers to 
/13 cm, assuming 30% thermal contribution 
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5. Discussion 

5.1. Cosmic-ray transport 



As found in Sect. l4.5l the helical magnetic field lines are con- 
fined to the walls that border the outflow cones, which appear as 
filaments when observed in projection. Therefore, we can study 
the cosmic-ray transport both, parallel to the filament direction 
and also perpendicular to it. 

Galactic cosmic-ray electrons with energies of a few GeV 
lose their energy mainly by synchrotron and inverse Compton 
(IC) radiation, and by adiabatic losses. Other losses like ioniza- 
tion ones can be neglected for the typical conditions of the in- 
terstellar medium. The conditions inside the n uclear region are 
differ ent, and ionization losses may dominate dThompson et al.l 
I2006I) . but we only investigate relativistic electrons here that 
have escaped from the nuclear region. 

The combined synchrotron and IC losses for a single cosmic- 
ray electron is given by 



dE 4 

d7 = ~r TC 



m e c- 



(t/ph + U B ), 



(3) 



where f/ p h is the photon energy density, Ub the magnetic field 
energy density, cr x = 6.65 x 10" 25 cm 2 is the Thomson cross 
section, and m e = 511 keV c~ 2 the electron rest mass. The time 
dependence of the energy is E(t) = Eq/(1 + t/t e ), where the 
electron lifetime is defined as the time at which the electron has 
lost half its initial energy 
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erg cm 



yr. 



(4) 



Assuming that the cosmic-ray electrons emit most of their en- 
ergy in synchrotron emission at the critical frequency, the elec- 
tron energy can be expressed by: 



B 



G)) 



1/2 



GeV. 



(5) 



I0 ergcirr 3 for B = 50 and U ph = 4.6 x 



To calculate the electron lifetime in the filaments we use 
U B = 1.0 x 10 
1 1 erg cm 

densities of the nucleus at a height of 300 pc using cylindrical 
symmetry. The IC and synchrotron losses are therefore compa- 
rable. The electron lifetimes are 2.5 (/120 cm), 1.4 (/16 cm), and 
1.1 x 10 6 yr 013 cm). 

For dominating radiation losses, the cosmic-ray electron 
scaleheight is a function of the frequency 



(/ -(l-a)/4 

h e oc y~ ' 5 (convection), 



h e oc v - u ~ a ' /4 (diffusion) 

-0.5 



(6) 



depending on whether we have diffusion or convection as the 
transport mechanism. These scaleheights are derived as follows. 
We note that in the equipartition case we have h e = 2/z syn (Eq. 1 1 
in Paper I). For diffusion, a is the energy dependency of the diffu- 
sion coefficient, usually assumed to b e a « 0.5 for Kolgomorov 
type turbulence dSchlickeiserl 120021) . For the case of convec- 
tion, the scaleheight is proportional to the electron lifetime with 
f e oc y~ 1/2 . In the case of diffusion we have /i 2 /f e = k, and with 
k oc E a and t e oc l/E we get the above relation. 

The observed frequency dependence of the scaleheight in the 
filaments is h e oc y 0±016 ; which is not consistent with convection 
in a galactic wind. Diffusion along the filaments would require 
a — 1 + 0.64, which can not be excluded a priori. However, 



diffusion should result in a Gaussian intensity profile, while the 
observed one is clearly exponential. Furthermore, the filament 
width of « 300 pc (Sect. l3.4l i is larger than its scaleheight. If 
cosmic -ray transport is responsible for the shape of the filaments, 
we would expect that the filament scaleheight is greater than the 
filament width, because diffusion along the field lines is much 
faster than perpendicular to them. We conclude that the cosmic- 
ray electrons do not propagate along the filaments, if radiation 
losses dominate. We give an alternative explanation in Sect. 15. 31 

The filaments have a Gaussian profile, as expected for diffu- 
sion, and their widths were determined in Sect. 13. 41 as a function 
of frequency as oc v -° 13 ± 018 This is consistent with diffusion 
perpendicular to the filaments with a = 0.5 ± 0.7, close to the 
expected behaviour, albeit with a large error interval. It seems 
therefore straightforward to assume that the cosmic -ray elec- 
trons of the surrounding medium diffuse from outside into the 
filaments. This explains why the filaments have a scaleheight 
nearly independent of the frequency, but their width decreases 
with increasing frequency (but see also Sect. 15.31 ). Summarizing, 
we can explain the observations by magnetic filaments with a 
width of about 40 pc (or less) at the boundary of both outflow 
cones. The perpendicular diffusion coefficient across the fila- 
ments is k ± = 1.5 x 10 28 cm 2 s- 1 ■ £(GeV)°- 5±0 - 7 . 

We can compare this result wi t h theo retical expectations as 
the recent paper by IShalchi et al.l d2010l) . who used numerical 
calculations to get an estimate for the perpendicular diffusion 
coefficient in a turbulent interstellar medium. Our cosmic-ray 
electron energy is a few GeV, so that for protons the magnetic 
rigidity R - (p ■ c)/ (Z ■ e) is R « 5 x 10 12 Volt. For this rigidity 
IShalchi et alJll2010h give a perpendicular diffusion coefficient of 



K x = 3 x 10 cm s , a factor of 10 smaller than what we find. 
However, the theoretically expected parallel diffusion coefficient 
is a factor of 100 larger. A small amount of turbulence in the 
filaments is able to explain our observations. 

5.2. Energetics in the nuclear outflow 

The nuclear outflow in NGC 253 appears to be a powerful ver- 
sion of the outflow in our own Milky Wa y. The supernova rate in 
NGC 253 lies between 0.03 and 0.3 yr 1 dUlvestad & Antonuccil 
1997; L enc &Tingavll2006l) and is thus much l arger than for 
the Ga lactic centre with (0.02 - 0.08) x 10~ 2 yr" 1 dCrocker et al.l 
2011a). Assuming that the outflow is driven by the starburst, the 
energy input in NGC 253 is thus a factor of 100 higher than in 
the Milky Way. Nevertheless, the Milky Way is able to launch 
a powerful galactic wind from its central region, as indicated by 
the strong d eficiency in the radio emission compared to the FIR 
luminosity dCrocker et al.l 1201 lbl) . Also, the gamma ray emis- 
sion is deficient compared to the star formation and supernova 
rates, which can be explained by the advection of cosmic rays in 
the wind, escaping into the halo . 

iDobler & Finkbeinerld2008l) have detected a microwave haze 
in WMAP observations towards the Galactic centre, which 
they interpreted as non-thermal synchrotron emission from 
cosmic -ray electrons. FERMI-LAlT observations have revealed 
two giant gamma-ray bubbles extending to a Galactic lat- 
itude of 50°, and containing an energy of 10 54 - 10 55 erg 
dDobler et aUl2010HSu et afll2010h . These discoveries were pre- 
ceded by the finding of a bipolar outflow in X-ray emission by 
iBland-Hawthorn & Cohenl d2003l) . estimating a similar amount 
of energy released. To create these bubbles, the current release 
of energy by star formation is not enough. However, the bub- 
bles have an age of 10 7 yr, so that a starburst in the last lOMyr 
could explain them. The alternative possibilities are past accre- 
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tion events of tidal debris of disrupted stars by the central black 
hole , which can a lso explain the presence of hot X-ray emitting 
gas (iTotanill2006[) . 

In the Galactic centre there are lobes of radio conti nuum 
emiss ion with a diameter of 1 lOpc and a height of 165 pc dLawl 
1201(H) . The Galactic centre lobes look very similar to the struc- 
ture we observe in NGC 253, as can be seen by comparing 
Fig. 1 in ILawl (12010b with our Fig.|4] although the filaments in 
NGC 253 are a factor of 2 larger. The equipartition field strengths 
between 40 and 100 //G are similar between the two galaxies. 
We note that we can trace the extension of the filaments in 
NGC 253 to a much larger height of 600 pc (Pig.[7j, likel y be- 
cause our beam size of 150pc is much larger than that of lLawl 
(l2010l) . Observations of the Faraday rotation favour a strong ver- 
tical m agnetic field in the Galactic centre lobe region (ILaw et al.l 
l20Tlh . 

We conclude that the radio continuum morphology of the 
Galactic centre and the nuclear starburst in NGC 253 are simi- 
lar, despite their current difference in star formation rate. As the 
nuclear outflow in NGC 253 can be explained by the starburst 
model alone, the similarity may hint that the Milky Way bipolar 
outflow seen in X-ray and gamma-ray emission is also generated 
by a past starburst, rather than the past activity of the Galactic 
central black hole. Gamma rays are now also detected from the 
nucleu s in NGC 253 by H.E.S.S., a Cerenkov air shower tele- 
scope (lAcero et al .120091) . According to their results, the nucleus 
is not a good "calorimeter", indicating a loss of cosmic rays, pos- 
sibly enlightening a bipolar outflow seen in gamma rays as in the 
Milky Way. 

5.3. Radio-FIR correlation 

As derived from IRAS flux densities, the nucleus of NGC 253 
has a FIR luminosity of FIR/(3.75 x 10 12 Hz) = 4.9 x 
10~ 8 erg s" 1 cirT 2 . Converting this into the radio (spectral) lu- 
minosity at 1.4 GH z, using the radio-FIR correlation (e.g., 
ICondon et ai]|2002l) . we find S 1.4GHz = 6.5 Jy. Our observed ra- 
dio luminosity of the nucleus of 2 Jy (Paper I) thus falls a factor 
of 3 short of the expectation. This impli es that the nucleus is 
not a "calorimeter" in the radio regime <IVoUdll989h . and that 
non-synchrotron losses of the cosmic-ray electrons could be im- 
portant. 

Such a process, which we expect to be strong in a galactic 
wind, is adiabatic cooling in an accelerated flow. If cosmic rays 
are advectively transported, the adiabatic loss time is 



4d = 3 (S 



(7) 



where dv/dz is the velocity gradient in the vertical direction. 
We note that the adiabatic loss time does not depend on the 
observing frequency, in contrast to synchrotron and IC losses. 
If adiabatic losses are dominating, we could naturally explain 
the independence of the synchrotron scaleheight from the ob- 
serving frequency (Sect. l3.4l i. As an upper limit for f ac j we take 
the cosmic-ray electron lifetime at A3 cm of 1 Myr, as derived 
in Sect. 15. II We use a cosmic -ray bulk speed of 300kms _1 as 
derived fr om the halo scaleheights (P aper I), which was con- 
firmed by IZirakashvili & Volkl d2006l) for the nuclear outflow. 
Assuming linear acceleration (with respect to z), we find that the 
cosmic rays are accelerated between and lOOpc height. 

We can also estimate the electron scaleheight as h e = v • f a d 
in the advection flow. The derived value of h e = 300 pc fits 
to the synchrotron scaleheight of 150pc, if the magnetic field 



scaleheight is much larger, such as the 600 pc expected for en- 
ergy equipartition. Since our estimated electron scaleheight fits 
the observations, the electron lifetime of 1 Myr is realistic. In 
this case, the adiabatic losses would be comparable to the com- 
bined synchrotron and IC losses. As the last two have about the 
same magnitude (Sect. 15. il l, we can explain a radio deficiency 
of a factor of 3. We conclude that adiabatic losses in the nu- 
clear outflow can explain the dim radio luminosity of the nu- 
cleus. But other processes may still contribute, such as ioniza- 
tion and bremsstra hlung losses of the cosmic rays in the nucleus 
(lLacki et alJl201lh . 



5.4. Collimating the nuclear outflow 

The SE outflow cone has an opening angle of 26° as Ha and 
X-ray data show (IStrickland et alJl2000t Bauer et al.ll2007l) . The 
X-ray data from the NW cone suggest a similarly small value al- 
though with a larger uncertainty. Interaction with the surround- 
ing interstellar medium from the galactic disc can explain such 
a small cone opening angle. Figures [3] and [4] show that the fila- 
ments of the radio emission border the outflow cones visible in 
the hot X-ray emitting gas. The Ha emitting gas borders again 
the X-ray outflow. In Fig. [17] we present a multi-wavelength 
composite of the central region of NGC 253, which nicely shows 
the interplay between the different components of the interstellar 
medium. The radio filaments clearly surround the outflow cones 
at least near the base. Therefore, it is useful to ask whether the 
magnetic field can coiiimate the outflow cone? 

To an swer this question w e make an estimate of the energy 
densities. Bauer et a D (l2007l) obtained an electron density of 
« e = 0.025 cm~ 3 at a temperature of 0.38 keV in the SE cone, re- 
sulting in an energy density of 1 .5 x 1 0" 1 1 erg cm~ 3 . The N W cone 
is brighter with a temperature of 0.54 keV so that here the en- 
ergy density is 3.0 x 10~ n erg crrT 3 (assuming also an increased 
electron density). The minimum estimate of the magnetic field 
strength in the filaments is 50 fiG resulting in an energy density 
of 1.0 x 10~ 10 erg crrT 3 . This is higher than the thermal energy 
density, so that the magnetic field can explain the collimation of 
both outflow cones. 



5.5. Galactic dynamo 

The parity of the magnetic fields can be used to test models of the 
galactic dynamo. The expectation is that a mean-field dynamo 
operating in a thin disc gives rise to a quadrupolar, even-parity 
magnetic field, where the direction of the azimuthal magnetic 
field is the same above and below the plane, but the direction 
of the vertical mag netic field componen t reverses with respect 
to th e plane (e.g., iKrause etail Il989bh lHaverkorn & Heesenl 
1201 ll) . This al so applies to a thin disc e mbedded in a quasi- 
spherical halo dMoss & Brandenburg! 19921) . However, a dynamo 
operating in a quasi-spherical halo should give dipolar, odd- 
parity magnetic fields, with an azimuthal magnetic field with 
reversing direction across the Galactic plane, while the vertical 
field is directed in the same sense above and below the plane 
dSokoloff & Shukurovll 19901) . So far the observational evidence 
prefers even fie ld parities in galaxie s, based on the RM patte rns 
in the disc (e.g jKrause et al.ll989allbL|Tabatabaei et al.ll2008l) or 
the asymmetry of polarized emission (IBraun et al ■1 12010b . These 
RM measurements of moderately inclined galaxies are, however, 
no stringent clue, because the RM structures of an even and odd 
parity in the disc fields are very similar in our synthetic polariza- 
tion and RM maps (Figs. A. 1 and A. 2 in Paper II). 
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Fig. 17. Three-colour comp osite, multi-wavelength view of the central region in NGC 253. Red, green, and blue indicate Ha from 
IWestmoquette et all d201 ll) . A20 cm radio continuum, and Chandra soft X-ray (M. Hardcastle, priv. com.), respectively. 



RM measurements in halo magnetic fields, which allow bet- 
ter distinction, are more difficult due to sensitivity and overlap 
with the disc field. Suc h measurements d o also favour quadrupo- 
lar, even-parity fields (ISoida et al.ll201 ll) . In Paper II we argued 
that NGC 253 is another candidate of an even disc field, albeit 
limited by the low resolution. This paper gives much better ev- 
idence for an even field in the disc beyond about 2kpc radius. 
High spatial resolution is crucial to distinguish between even 
and odd parity of the disc magnetic field. The synchrotron scale- 
height of the thin disc is about 300 pc (Paper I), so that we need 
a comparable resolution, like 6" (D = lOMpc), to resolve an 
RM reversal across the midplane. Such high angular resolutions 
with sufficient S/N in polarization now become feasible with the 
advent of new broadband correlators at various radio interfer- 
ometers. We expect that future observations will be better able 
to detect odd disc fields in galaxies, some of which may have so 
far been classified as even. 



sity within a radius of as 2kpc. Anisotropic magnetic fields 
can be generated from random, isotropic fields by shear and 
compression. M31 is a galaxy that is dominated by a co- 
here nt axisymmetric magnet ic field field from its RM struc- 
ture dBerkhuiisen et al.l 12003). In contrast, M51, a grand-spiral 
galaxy with a prominent spiral magnetic field , shows almost no 
large-scale RM pattern (Fle tcher et al.ll20TTl) . This galaxy is a 
prime example of an anisotropic disc magnetic field. The late- 
type spiral galaxy NGC 6946 has both regular magnetic fields 
between the spiral arms, in the so-called "magnetic arms", and 
anisotropic magnetic fields (Beck 2007). Barred galaxies, such 
as NGC 1097 and NGC 1365, wi th their strong sh earing flows 
host dominant anisotropic fields (Be ck et aT1l2005l) . In light of 
these results, it comes as no surprise that the ordered disc field 
in NGC 253 is anisotropic within « 2kpc radius, where its bar 
resides (Sect. l3.U . 



We detected anisotropic magnetic fields as the dominant con- 
tribution to the ordered disc field observed in polarized inten- 
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5.6. Helical magnetic field in the outflow cone 

The helical field can be explained with help of trajectories of par- 
ticles in the outflow. The conductivity of the interstellar medium 
is very high, so that the magnetic field lines are essentially frozen 
into the plasma. The field lines are taken with the gas in the 
outflow similar to the solar wind dParkerlll958t IWeber & Davisl 
The magnetic field is anchored in the disc between a ra- 
dius of about 130 - 170pc (Fig.lT5T>. The cone opening angle 
of 26° (45° at larger heights) means that the field is transported 
to larger radii farther away from the disc. If we assume that no 
significant interaction with the surrounding medium occurs, the 
angular momentum of the plasma in the outflow is conserved. 
From the perspective of an outflowing particle, the anchor point 
rotates faster than the particle. The trajectories of the particles 
can be described by a spiral, similar to water drops in a garden 
sprinkler. The spiral orientation is given by the sense of rotation 
and is identical to that of the spiral arms, i.e. trailing arms. The 
vertical component of the magnetic field has to point away from 
the disc in order to explain the observed RM structure, which 
agrees with our results in Paper II for the SE halo. 

We can compare the actual shape of our helix with this 
model. The rotation curve of the galaxy is rising linearly out 
to a radius of 300 pc, where it tailo rs off to a rotation speed of 
about 200kms _1 (ISorai et al.ll2000l) . The rotation speed at the 
anchor point of the filaments is thus « 100 km s , resulting in a 
rotation period of 0.9 x 10 7 yr. We measured that the azimuthal 
component is about equal to the vertical component at a height of 
1200 pc (Sect. l4.5l l. The circumference at 150pc radius is 940 pc, 
similar to the height, so that we can expect an equal strength of 
the vertical and azimuthal components. During one rotation pe- 
riod, the particles in the outflow should thus be transported to a 
height of 120 pc. The outflow speed in the nuclear outflow is 
« 300 km s _1 ( Wes tmoquette et al.l201 ll) at a reference height of 
1 kpc, so that assuming constant acceleration the average speed 
is 150kms~'. To reach a height of 1200 pc a particle would take 
a time of 0.8 x 10 7 yr, which is in good agreement with the rota- 
tion period. 

Our simple model of the magnetic field lines frozen into the 
outflow plasma can therefore explain the observed winding-up 
of the field lines into a helix. We also note that the Alfven speed 
Va = Bj -\J4-Kp is 30kms _I with n e = 2cirT 3 and B — 21 /yG. 
This is only 10% of the outflow speed at a reference height of 
1 kpc, so that the Alfvenic point would be with lOOpc close to 
the disc. The stiff field lines co-rotate with the disc within the 
Alfvenic point. This means that the field lines start almost im- 
mediately winding up in the outflow, which accords with our 
observations. 

Although the nuclear outflow in NGC 253 is thought to be 
driven entirel y by a starburst and no t by an active galactic nu- 
cleus (AGN, Brunthaler et al.l l2009t) . it is worthwhile briefly 
comparing the outflow structure with those of AGNs. There 
are a number of radio polarimetric studies of AGNs, employ- 
ing very long baseline interferometry and studying their mag- 
netic fields using the sy nchrotron emission of high ly energetic 
electrons in the jet (e.g jGabuzda et alj|1992ll2000h . They show 
that the magnetic field orientations are mainly transversal to the 
jet direction. This is thought to be caused by compression in 
shock waves, where the cosmic -ray electrons are accelerated, 
and thus illuminate the "knots" in the jets where the shocks oc- 
cur. However, the alternative explanation is that the transverse 
fields are generated by the diffe rential rotation in the acc retion 
disc around the black hole (e.g.. lContopoulos et alJl2009b . This 



would result in a helical magnetic field farther away from the 
accretion disc. 

One signature of a helical magnetic field in the jet of AGNs 
should be a reversal of the RM across the jet, as observed by 
us in the nuclear outflow of NGC 253. Studies of the RM in jets 
of AGNs are notoriously difficult, given that the jets are mostly 
only barely resolved and the thermal electron densities are low 
dGabuzda et al.ll2004t iKharb et alj|2009t |Q' Sullivan & Gabuzdal 
|2009t). However, these authors find a strict dependence of the 
polarization angle on A 2 , suggesting that the jets form a Faraday 
screen and that the observed RM gradients are genuine and not 
instrumental. If the view onto the jet is nearly "side-on" (90° to 
its axis) in its rest frame, the magnetic field would be transversal. 
If the jet is observed at a different angle to 90°, we should see a 
central region with a transversal magnetic field, whereas the field 
should be along the jet direction at the boundary. Such "spine + 
sheat" structures have been i ndeed observed (e.g jAttridge et all 
ll999t|Pushkarev et al.l2005h . further supporting the existence of 
helical fields in AGN jets. 

6. Conclusions 

We used polarimetric radio continuum observations to study the 
nuclear outflow in NGC 253. Four prominent filaments were de- 
tected, which we interpreted as the walls of the two nuclear out- 
flow cones seen in projection. The boundary of the NW outflow 
cone is not seen in Ha and X-ray emission, probably due to 
heavy absorption. These are our main conclusions 

1 . The filaments have a scaleheight of 150 + 20 pc, independent 
of wavelength. The cosmic -ray electron lifetime changes by 
a factor of 2.3 between AA 20 and 3 cm, so that the scale- 
height should change if the electrons are subjected to domi- 
nant radiation losses. Dominating adiabatic losses could ex- 
plain this behaviour, where cosmic rays and magnetic fields 
are advected together in the accelerating nuclear outflow. 

2. The filaments have a smaller perpendicular extent for smaller 
electron lifetimes, so that perpendicular electron diffusion 
plays a role. This suggests that we are looking at filamen- 
tary magnetic fields. The perpendicular diffusion coefficient 
is/c x = 1.5 X 10 28 cm 2 s- 1 • £(GeV) a5±0 ' 7 . 

3. The filaments have a total magnetic field strength of B tot = 
46 + 10 iiG and ordered field strength of B old = 21+5 j«G 
The field strengths could be higher, if the cosmic rays are not 
in equipartition with the magnetic field. 

4. The magnetic field strength in the filaments is high enough to 
counteract the pressure in the nuclear outflow. The magnetic 
field can thus collimate the outflow and explain the small 
opening angle of « 26° by interaction with the interstellar 
medium. The magnetic field is aligned along the filaments, as 
expected for the case of compression due to a lateral pressure 
gradient. 

5. In the SE outflow cone, the magnetic field points away from 
the disc in the form of a helix. The magnetic field is frozen 
into the plasma of the nuclear outflow, so that the angular 
momentum is conserved and the field wound up. 

6. The ordered magnetic field in the disc observed in polarized 
intensity is anisotropic within a radius of about 2 kpc. At 
larger radii a regular field with an even axisymmetric pattern 
is found, as already deduced in Paper II. 
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Appendix A: Synthetic polarization maps 

In lHeesen et al ] (l2009bl) we used a conical halo field with a filled 
cone while in this paper we use a conical field that is only differ- 
ent from zero inside the walls of the outflow as shown in Fig. [15] 
The following equations are therefore identical for both cases. 
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We stress that, even though we denote the vertical magnetic field 
component with the superscript "halo", this field is not a general 
large-scale halo field but is only concentrated in the filaments. 
We do not use any turbulent magnetic field, but instead restrict 
ourselves to ordered magnetic fields. 

We integrate the source structure in a three-dimensional vol- 
ume with Cartesian coordinates (x, y, z), where the x-axis is par- 
allel to R.A., y is parallel to the line-of-sight, and z parallel 
to Dec. We also define the coordinates (x',y',z') that are the 
Cartesian coordinates with respect to the galaxy, x' is parallel 
to the major axis, y' is parallel to the minor axis, and zf parallel 
to the rotation axis. We can relate both coordinate systems by 
linear transformations first rotating around the inclination angle 
; and then rotating it around the position angle. For convenience 
we use i = 90° - ;, which is the angle between y and y' , and 
£ = 90° - p. a. which is the angle between z and z' ■ The ref- 
erence coordinates (x re f , y re f , z re f) are the central coordinates of 
the galaxy. Moreover, we use cylindrical coordinates for the ob- 
server (p, <p, z) and in the galaxy (p', zf). 

The frame of reference Cartesian coordinates of the galaxy 
can be expressed in the observer's as 

x = cos(£) • (x - x re f) - sin(£) • (z - z ref ) (A. 1) 

y' = - sin(£) • sin(t) • (x - x ief ) + cos(t) • (y - y ref ) 

- sin(i) • cos(£) • (z - Zref) (A.2) 

z = sin(£) ■ cos(i) ■ (x - x ref ) + sin(i) • (y - y K{ ) + 

COS(t) • COS(£) ■ (Z - Zref)- (A.3) 

For the strength of the magnetic field in the disc we use a 
Gaussian profile with respect to the galactocentric radius and 
an exponential vertical decrease with a scaleheight /is.disc- For 
NGC 253 we used a two-component Gaussian distribution as a 
function of p', which is not shown here for simplicity. The pa- 
rameter cr can be expressed by the full width half mean (FWHM) 
as cr = FWHM/(2 Jl log(2)) in this example and all following 
cases where we use a Gaussian distribution: 



gdisc _ ^disc 



exp 



exp 



-P 



- A 



^^"disc ^ 



(A.4) 



where we take the absolute value |z'| of the vertical distance to 
the galactic midplane. The magnetic field strength in the halo 
(filament) can be expressed as the disc field by 



5" 



exp 



• exp 



^^"halo 



(A.5) 



are the strengths of the ordered field in the 
disc and halo. The z'-component of the halo (filament) field is 



b: 



^halo 



z,halo 



Vl.0 + tan 2 (ff) + tan 2 (j8) 



(A.6) 



where for even parity "+" is used if z' < and for odd parity "+' : 
is used for all z'. We can now derive the ^'-component as 



^V.halo _ ± ^c,halo ' 



tan(of), 



(A.7) 



where "+" is used if z' < 0. The angle a gives the azimuthal 
component of the halo (filament) field and B is the angle with 
respect to the rotation axis, i.e. half of the cone opening angle. 
Now 



^haio = ±#ha,o-tan(/?), 



(A.8) 



where "+" is used if z' < 0. Now 



B' 



^disc 



0',disc 



1.0 + tan 2 O) 



(A.9) 



where for even parity "— " is used for all z' and for odd parity 
"+" is used if z' < 0. The azimuthal disc component is 



•Bp',di sc - +tan(y) ■ -00- .disc, 



(A. 10) 



where - is used for an inwards pointing field. For the following 
we define <p as 



(p = atan2(y', x'), 



(A.ll) 



where the numerical function "atan2" is equivalent to atan, but 
returns the angle in all four quadrants (360°) rather than only 
between -90° and 90°. Transforming back from the cylindrical 
coordinate system of the galaxy into the Cartesian coordinate 
system of the observer we first need to transform to the Cartesian 
system in the galaxy: 



(flf c )' = (fl* sc )' • cos(^) - (B** c )' ■ sin(^) 

(sf °y = (sf y ■ sin(^) + (Bfy ■ co^y 



(A. 12) 
(A.13) 



Transforming above equations into the Cartesian system of the 
observer, 



Bf c = cos(£) • (Bf c )' - sin(£) • sin(i) ■ (Bf c )' + 
sin(£) • cos(i) ■ (Bf sc )' 



Bf c = cos(t) • (Bf c )' + sin(i) • (Bf sc )' 



Bf sc = - sm(0 ■ (Bf c )' - sin(t) • cos(£) • (Bf c )' + 
cos( t ) ■ cos(£) ■ (Bf c )'. 



(A. 14) 
(A. 15) 

(A. 16) 



The same equation applies to the halo (filament) magnetic field 
that is not shown here. Now we have the magnetic field compo- 
nents in the Cartesian system of the observer. We define the local 
polarization angles for the disc and halo field as 



<A disc = atan2 {-(Bf c )', (B dlsc )'} - 90° 
<A hal0 = atan2 {-(Bf lo )', (Bf lo )'l - 90°. 



(A. 17) 
(A. 18) 



The 90° correction is to transform the polarization angle into 
the magnetic field orientation. The Faraday rotation shifts the 
polarization angle if we go a step Ay along the line-of-sight: 



A^(y + Ay) = Ai[r(y) - 0.81 ■ A ■ n e ■ exp 

I n xdisc 

P -(Bf c )' + 



(A. 19) 



The "-" is because we are going in positive y-direction, whereas 
the polarized wave is going from the source to the observer. The 
ratio of the regular to the ordered field describes the anisotropy 
of the magnetic field, with B reg < B 01 d- 
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(a) 




(b) 




Fig. A.l. Coordinates used for the synthetic polarization maps, (a) View along the line-of-sight in the direction of the y-axis. (b) 
View from the side in the opposite direction to the jc'-axis, where the observer is located at the left. 



We can now also calculate the Faraday screen RM if the Q(x,z) and U(x,z) are written out as a fits file for which we 



galaxy would be a screen: 



used cfitsi 



?M screen (y + Ay) = RM scleen (y) - 0.81 ■ « e ■ exp| -j— 



exp 



(Bf sc )'+ 



i halo 



(Bfy \ Ay. 



(A.20) 



For the observed polarization angle we obtain with Faraday ro- 
tation included: 



(A.21) 



The cosmic-ray energy density «cr is scaled according to the 
equipartition condition as 



"CR - »CR,0 



' (Bf sc ) 2 + (Bf sc ) 2 + (Bf sc ) 2 
, « sc ) 2 + (-O 2 



(A.22) 



(B$ sc ) 2 + (B% dio ) 2 

The Stokes Q integration along the line-of-sight is 

Q,U(y + Ay) = Q,U(y) + n CR - 

{((Bf^ 2 + (Bf-) 2 )-co S (2^f dy ) 

+ ((Bf °) 2 + (B? al °) 2 ) ■ cos(2<:;: day )} , (A.23) 

6 CFITSIO is a library of C and Fortran sub-routines for read- 
where for Stokes U, "sin" is used instead of "cos". Stokes Q and ing and writing data files in FITS (Flexible Image Transport System) 
U have to be calculated at at least two frequencies in order to d ata format. CFITSIO is free software an d can be obtained at 
calculate the rotation measure. The two-dimensional arrays of http://heasarc.gsfc.nasa.gov/fitsio 



